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Abstract: 

Recently, Hellerman and Schnabl considered the dynamics of unstable D-branes in the 
, background of a linear dilaton. Remarkably, they were able to construct light-like tachyon 

solutions which interpolate smoothly between the perturbative and nonperturbative vacua, 
without undergoing the wild oscillations that plague time-like solutions. In their analysis. 



5h ' however, the full structure of the initial value problem for the nonlocal dynamical equations 

was not considered. In this paper, therefore, we reexamine the nonlinear dynamics of light- 
like tachyon condensation using a combination of numerical and analytical techniques. 
We find that for the p-adic string the monotonic behaviour obtained previously relied 
on a special choice of initial conditions near the unstable maximum. For generic initial 
conditions the wild oscillations come back to haunt us. Interestingly, we find an "island of 
stability" in initial condition space that leads to sensible evolution at late times. For the 
string field theory case, on the other hand, we find that the evolution is completely stable 
for generic choices of initial data. This provides an explicit example of a string theoretic 
system that admits infinitely many initial data but is nevertheless nonperturbatively stable. 
Qualitatively similar dynamics are obtained in nonlocal cosmologies where the Hubble 
damping plays a role very analogous to the dilaton gradient. 
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1. Introduction 

Nonlinear theories with infinitely many derivatives have come to play an increasingly im- 
portant role in theoretical physics; attracting interest both from string theorists and cos- 
mologists. The string theory application of such equations that has stimulated the most 
interest is that of understanding the dynamics of tachyon condensation in string field the- 
ory (SFT) Ijl], H and also in related toy models, such as the p-adic string Q. The dynamical 
process of the tachyon field rolling from the unstable maximum of its potential to the true 
vacuum is expected to give a time-dependent description of the decay of unstable D-brane 
configurations in string theory (see and references therein). A long-standing puzzle has 
been the observation that generic tachyon solutions in flat space do not roll monotonically 
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from the unstable maximum of the potential to the true vacuum of the theory, as would be 
expected on physical grounds. Rather, the tachyon field undergoes wild oscillations at late 
times ||5[ ; a manifestation of the Ostrogradski instability |p that plagues higher derivative 
theories (see Q and Q for a more modern discussion). The string theoretic interpretation 
of this peculiar behaviour is a subtle problem and there are a number of proposed expla- 
nations in the literature (discussed in more detail below). (See Q for further discussion of 
tachyon dynamics in SFT and related theories and see |jl^, |lT[ for significant progress in 
understanding the vacuum structure of SFT.) 

The wild oscillations of the tachyon field arise because the higher derivative structure 
of SFT allows for extra (more than two) initial data in the solutions of the field equations. 
These extra data can be interpreted as a tower of new physical states, in addition to the 
usual tachyonic excitation, which contribute to the kinetic energy with indefinite sign |13|. 
In the case of SFT and the p-adic string, these spurious extra states have complex mass- 
squared and behave as an admixture of ghost and non-ghost field (we will refer to such 
excitations as ghost-like or "quintom" in the text). The presence of extra ghost (or ghost- 
like) degrees of freedom is a rather generic problem for higher derivative theories; see Q 
for a detailed discussion. In any theory with negative kinetic energy the Hamiltonian is 
unbounded from below and the system will become arbitrarily excited at late times for 
generic choices of initial data. The kind of peculiar time dependence obtained for rolling 
tachyon solutions in flat space is quite typical of higher derivative instabilities. 

Recently, Hellerman and Schnabl |14] made significant progress in studying the dy- 
namics of brane decay by explicitly constructing tachyon solutions that roll smoothly from 
the perturbative vacuum to the true vacuum, as a function of light-cone time, x"*". This 
progress relied on turning on a linear dilaton background that violates energy conservation 
and provides a source of friction for the tachyon dynamics. Hellerman and Schnabl con- 
sidered light-like tachyon profiles in SFT, p-adic string theory and also vacuum string field 
theory (VSFT), finding similar dynamics in all three cases. Despite their significant result, 
however, the analysis of does not explicitly consider the role of initial conditions near 
the unstable maximum. Without a complete understanding of the initial value problem 
it is impossible to assess the stability of such solutions. Here we point out that the so- 
lutions presented in |14] result from the choice of one particular kind of initial condition 
from an infinite set of possibilities. Furthermore, we will show that for the cases of the p- 
adic string and VSFT the monotonic behaviour obtained by Hellerman and Schnabl is not 
generic in the sense that it relies on the special choice of initial conditions made. For more 
general choices of initial conditions the wild oscillations associated with the Ostrogradski 
instability come back to haunt us. On the other hand, we will show that the equation 
obtained in [0] for SFT in a level zero truncation leads to completely stable dynamics, 
even when the extra initial data are taken into account! This behaviour is confirmed by 
fully nonlinear numerical simulations, and provides a remarkable example of an interacting 
nonlocal theory that admits infinitely many initial data but whose evolution is nevertheless 
completely stable. These surprising dynamics arise because the friction coming from the 
dilaton gradient efficiently damps out the oscillations of the tachyon field. 

A final surprise awaits us. Further studying the p-adic and VSFT case, we discover that 
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despite the instability associated with generic initial conditions, the set of initial conditions 
leading to stable evolution is not of measure zero. This implies that even in these cases 
there is an "island of stability" in the initial condition space which leads to a well behaved 
evolution. 

Nonlocal theories motivated by SFT have recently attracted interest also from cosmol- 



ogists [15]-[36| due to a wide array of novel cosmological behaviours. Of particular interest 



are recent efforts to construct inflationary solutions in nonlocal theories [^]-|36]. The first 
realisation of nonlocal inflation was in the context of p-adic inflation [^] . This model can 
support inflation even when the potential is naively very steep, a behaviour that was found 



to be a rather generic feature of nonlocal inflation in pH] (see also [p3[) and was verified 



using numerical analysis in |36|. Moreover, nonlocal inflation is also one of the rare models 
that predicts large non-gaussianity in the cosmic microwave background ||3^, Here we 
point out that the nonlinear dynamics of nonlocal inflation are strikingly similar to the 
light-like tachyon models considered by Hellerman and Schnabl. Hence we expect that our 
analysis should have implications also for the stability of nonlocal cosmologies. 

During the course of our investigation we uncover a number of results of general interest 
to the study of nonlinear theories with infinitely many derivatives. For instance, we will 
argue that the recipe of mixing friction and constraints on the initial conditions provides 
a very generic prescription for constructing stable solutions in infinite order theories. We 
believe that this work should be useful in guiding future searches for stable solutions in SFT 
(and similar theories). Moreover, our analysis clarifies a number of issues concerning the 



mathematical structure of nonlinear infinite order differential equations. In [12| the initial 
value problem for linear constant coefficient equations with infinitely many derivatives was 
studied and a formalism was developed to exhaustively count initial data.^ These results 
were generalised to the case of variable coefficient equations (such as those that arise 
when studying nonlocal cosmological perturbation theory) in p^. However, the nonlinear 



problem presents a number of questions which could not be addressed in [12, 37 1. This 
work represents progress towards a complete understanding of nonlinear equations with 
infinitely many derivatives. 

The organisation of this paper is as follows. In section |2| we study the nonlinear dy- 
namics of the light-like p-adic string tachyon in a linear dilaton background, presenting 
non-linear analytic solutions which fully take into account the freedom in fixing initial 
conditions. In section |^ we study the analogous dynamics at level zero truncation in SFT. 
In this case analytical solutions are not available so we must turn to a numerical analy- 
sis. In section Q we describe our formalism for solving infinite order differential equations 
numerically, before applying our numerical approach to study the light-like SFT tachyon 
with generic initial conditions in section ^. In section ^ we comment on the stability of 
nonlocal cosmologies. We briefly review some proposed stringy interpretations of the wild 
oscillations in rolling tachyon solutions in section 0. Finally, in section ^ we conclude. 



^See also different approaches to the initial value problem and for mathematical analysis of 

p-adic and string field equations. 
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2. Light-Like Tachyon Condensation in p-adic String Theory 



2.1 Set-Up and Equation of Motion 

We begin our investigation of the dynamics of light-like tachyon condensation by studying 
p-adic string theory Q coupled to a linear dilaton profile. We employ the action proposed 
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(with 



where p is a prime number that characterises the world-sheet coordinates, a' 
nis the string mass), Qp is related to the open string coupling and <j) is the (dimensionless) 
tachyon field. Although ( p.lD is not meant to be construed as a realistic model of string 
theory, we will spend some time studying this theory because it is analytically tractable 
and provides an excellent playground for studying the nonlinear dynamics of infinite order 
theories. Following [14| we will work in terms of light-cone coordinates = (x" it x^)/^/2 



so that the metric takes the form ds'^ = —{dx^Y + {dx'^Y + dy- dy = —2dx~^dx~ + dy ■ dy. 
The D — 2 transverse coordinates y will play no role in the ensuing analysis. The dilaton 
background is taken as 



= Vf.x^' = -V+x^ - V-X+ + V-y. 
For a light-like field <j) = (j){x~^) one obtains the equation of motion 



X + a V Inp) = p 



(2.2) 



(2.3) 



Note that this finite difference equation can trivially be re-written as a pseudo-differential 
equation 

p'^'^^^+</.(x+) = p"'^'/2^P(x+) . (2.4) 

This equation admits constant solutions (f) = ^""'^^/[^(p^i)] and (j) = 0. The former corre- 
sponds to the unstable maximum of the potential (physically the state with a space-filling 
brane) while the latter is the true minimum (physically the nonperturbative vacuum with 
no D-brane). 

The action ( |2.lD was motivated by the striking similarity Eq. (^|^) bears to the equation 
of motion one derives for the tachyon in VSFT (with light-like ansatz and linear dilaton 



background). Equation (2.4) is also identical (up to a re-scaling of the field and space-time 
coordinates) to the friction dominated equation for the inflaton dynamics in p-adic inflation 
j32U.'^ Due to these similarities we expect our analysis of (2.4) to have relevance to VSFT 
and p-adic inflation also. 

2.2 Exact Analytic Solution 



Remarkably, equation (2^) admits an exact nonperturbative solution: 



ix 



(2.5) 



^Setting = and replacing — > t, 3Ho/2 one obtains exactly the equation that was studied 

in section 4.3 of IB2|. 
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where F(x~^) is an arbitrary smooth periodic function satisfying 

F{x+) = F{x+ + a'V+ Inp) . (2.6) 
It follows that F can be decomposed into a Fourier series as 




It is clear that the solution contains an infinite number of free parameters {a„,, 6„} which 
allow us to fix the state of the solution at some initial time (we will usually set x^ = 
in subsequent analysis). Hence an infinite number of initial conditions must be specified 
to give us a unique solution. As we will see, the first term in ( p.7D (proportional to ao) 
is associated with the tachyonic excitation while the the oscillatory terms (proportional to 
a„,6„ with n > 0) are associated with the presence of ghosts in the theory, since these 
modes contribute with indefinite sign to the total kinetic energy. We will refer to the n ^ 
states as ghost-like (or quintom) modes. 

Let us discuss the dynamics of the solution ( |2.5| ) in more detail. For any choices of an 
and 6„, the solution rolls from the unstable maximum in the asymptotic past x"*" — > — oo. 
For Oo > and a„ = 6„ = for n > the field rolls monotonically from the unstable 
maximum to the true vacuum of the theory. This solution corresponds to turning on only 
the tachyon in the initial state and in this case our solution ( p. 5]) matches the one derived 
in (also identical to the solution obtained in [^2| in a somewhat different context). 

On the other hand, considering oq = and taking any of the or bn to be non-zero 
leads to a wildly oscillating solution which never settles to the minimum, but rather the 
amplitude of the oscillations grow in an unbounded manner (similarly to the solutions of 
( |2.1| ) in the background of a constant dilaton This is not unexpected, since these 

modes of excitation are ghost-like and lead to vacuum instability. 

It is extremely interesting and indeed unexpected, however, to see what occurs if we 
choose to turn on some admixture of ghost excitation together with the well behaved 
tachyon in the initial state. Let us again take oq > 0. In this case we can choose the origin 
of time such that oq = 1 without loss of generality. Now let us arbitrarily take some of the 
an,bn to be non-zero, but sufficiently small in a sense which we will quantify shortly. In 
this case the solution again rolls away from the unstable maximum, but as it rolls down 
the potential the field undergoes some small oscillations, though it still settles down to the 
minimum at late times. Quantitatively, the condition 

F{x+) > (2.8) 

is sufficient to ensure that (/)(x^) — > as x"*" — > oo. Let us now consider increasing a„, 6„ so 
that ( |2.8D is violated. In this case the ghost excitations become dominant and (p oscillates 
wildly at late times. Note that these initial conditions are physically reasonable because (p is 
at the false vacuum in the asymptotic past x^ — > — oo. In Fig. I we illustrate the different 
behaviours of the solution (|2.5| ), which we have just discussed, for some representative 
choices of a™, 6^. 
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p-adic string 




Figure 1: The solution Eq. (|]^), as a function of x'^ . The curves are ai = (dashed hue), 

ai = 0.25 (soHd Une), Oi = 1.1 (dash-dotted Hne) with all other a„ and 6„ vanishing and ag — 1. 
For illustration we have set p — 3 and — 0. The x'^ axis is measured in units of a'V'^. 



Let us comment on the structure of the initial value problem for (^^), which is closely 
related to the (in)stability of the theory. In [12| a formalism was presented for exhaustively 
counting the initial data of linear, infinite order equations (see |37] for a generalisation to 
variable coefficient equations). One could use this formalism to i U-4 ) in ^ perturbative 



expansion around 



-a'Vy[2{p~l)] 



The solution so obtained reproduces the terms 



in a small- It expansion of (|2.5| ). Using the results of |12| then proves that the solution 
( |2.5| ) provides a complete solution of the initial value problem near the false vacuum. By 
continuity, we expect that this solution is complete also nonperturbatively. 



2.3 Behaviour of the Hamiltonian 



Using our exact solution (|2.5|) it is straightforward to explicitly construct the energy density 
as a function of (this is not simply a constant because the dilaton gradient violates time 
translation invariance). For simplicity we consider V~ = ^ = so that <I> = —V~^x~ and 
y2 = although we do not expect that our qualitative results depend on this restriction 
in any crucial way. The nontrivial components of the stress tensor are [0] 



'-++ 



a'lnp _<j, 



9l 



+ 



dCd+ [0(x+ + V+a'lnpC)] 



(2.9) 
(2.10) 
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Evaluating these on the solution (^]^) we have 



2 ^rr 



2 ^rp 



-2e-+/(«'^+)F(x+) 



-{p+l)e^ /(" )F(x+) 



2 p+1 

.2x+/{a'V+) 



(2.11) 



a'(y+)2 



exp 



j\^e-^e^*''-'''*'F(-^+^'V+^nz) [p {x+ + a'V+ z) + 



aV+F'(x+ + a Inz)] 



(2.12) 



where in evaluating T++ we have switched variables z = p'' in the integration. Remembering 
that T^jy transforms as a tensor the energy density is 



P 



(2.13) 

Since ( 2.11| ) could be written as T-| = —e^'^U{cj)) where U{(f)) = — it is natural to 

associate the second term in ( 2.13| ) with the potential energy and the first term with the 
kinetic energy. It can be readily verified that the kinetic term has indefinite sign and that 
the presence of some non-zero a„,6„ (with n> 1) in ( |2.7D leads to negative contributions 
to the kinetic energy. 

energy density in the p-adic string 
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Figure 2: The energy density g^e^ p as a function of a;"*" for the same values of p and a„, 6„ used 
m fig. 0. The a;+ axis is measured in units of a'V^ and, for illustration, we have set = 1 in 
these units. 

The behaviour of the energy density as a function of is plotted in Fig. |2|. In the 
case Qn = bn = for n > (the case considered by [jl^) we have monotonic behaviour of 
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p{x'^). Turning on the ghost-like n > modes the osciUations of (p lead to oscillations in 
p. If the contamination of ghost modes in the initial state is sufficiently large then these 
oscillations cause p{x~^) to cross zero and one would expect recollapse if gravity had been 
included. From the previous discussion of (l){x~^) one might expect that is a sufficient 
condition to keep the kinetic energy positive for all x^. However, detailed examination of 



(2.13) reveals that the actual condition is somewhat stronger: 

+ a'y+F'(2;+) > . (2.14) 

The fact that (|2.14| ) is a stronger constraint than ( |2.8D means that we can construct so- 
lutions where p{x^) crosses zero without prohibiting (l){x~^) from settling down to the 
minimum at late times. 

2.4 The Island of Stability 



The instability of the theory (2.4) is not surprising. This equation admits infinitely many 
initial conditions and hence Ostrogradski's theorem implies that the dynamics should be 
generically unstable. Of course it may be possible to construct well-behaved solutions by 
carefully choosing initial conditions, however, one would expect the set of initial conditions 
leading to non-pathological evolution to be a set of measure zero. This intuition is based 
on the observation that the Ostrogradski Hamiltonian is linear in the unstable directions. 
However, this intuition is incorrect in our case: the set of initial conditions leading to 
sensible evolution is not measure zero. Equation (|2.4| ) allows for a certain amount of 
ghost-like contamination in the initial conditions (quantified by taking a„,6„ to be non- 
zero) without prohibiting the tachyon from settling down to the vacuum at late times. This 
observation implies that there is an "island of stability" in initial condition space. This 
(perhaps) surprising behaviour arises because the dilaton friction can efficiently damp out 
the unstable growth of the ghost modes as long as they do not contaminate the initial state 
too strongly. 

It is worth emphasising that the island of stability is not a perturbative artifact. Our 
solution ( ^.51) is exact at the fully nonlinear level and, as we have argued, it provides a 
full solution of the initial value problem for the theory ( |2.4| ). It is tempting to argue that 
a theory like ( |2.4| ) could be phenomenologically viable once the initial data are suitably 
restricted. Such a restriction on the data could be elegantly implemented using the contour 
deformation prescription advocated in [12|.^ We would, however, like to emphasise two 



caveats. First, our analysis is purely classical. Quantum mechanically a solution which 
starts in the island of stability might be able to tunnel to an unstable configuration, perhaps 



though a negative tension instanton solution similar to those constructed in |44]. Second, 
we should emphasise that the stability of the solution ( |2.5| ) is closely tied to the light-like 
ansatz (j) = (j){x~^). On scales small compared to the dilaton gradient one can treat $ as a 



constant and (^.1]) should have all the usual instabilities |§| and wild oscillations 



^The expression "island of stability" is taken from ref. where similar behaviour was claimed for a 
finite order higher derivative theory. 



*The solution (2.5) provides an explicit example showing that this prescription can be applied at the 
fully nonlinear level. 
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It is interesting to compare our findings to the case of closed string tachyon condensa- 
tion. Purely light-like closed string tachyon solutions are stable against small perturbations 
of the initial conditions (in the sense that small perturbations of the initial conditions lead 
to qualitatively similar behaviour at late times) 46 1. On the other hand, [47| discussed 



the possibility of qualitative change of the solution under finite changes of the initial con- 
ditions. 

3. Light-Like Tachyon Condensation in String Field Theory 
3.1 Set-Up and Equation of Motion 

We now consider light-like tachyon condensation in the more realistic context of string field 



theory at level zero truncation. The action derived in [14| for the SFT tachyon in the linear 
dilaton background (given by Eq. ( |2.2| )) is 

(3.1) 



where go is the open string coupling constant and the constant \s K = 4/33/^. Imposing 
again a light-like ansatz on the tachyon field (j) = (j){x'^) we obtain the equation of motion 

(aV+a+ - 1) 0(x+) + -^0^(2;+ + 2a V+ \nK) = {). (3.2) 
This can be re-written in pseudo-differential form as 

(a'V+d+ - 1) K-2a'y+9+^(^+) ^ -^(j)'^{x+) = . (3.3) 

This equation admits constant solutions 1^ = and (f) = K^. The former corresponds to 
the unstable maximum while the latter is the true vacuum. 

Note that, as in the p-adic case, equation ( |3.3| ) is identical to what one would have 
for the time-like SFT tachyon in a de Sitter background in the limit of very large Hubble 
scale where one could take □ = -df - 3Hdt = -3Hdt. Hence, in some (limited) sense the 
dilaton gradient acts like an infinite source of Hubble friction. 

3.2 Perturbative Analysis 

Sadly, the equation ( p. 3D does not seem to admit a simple exact analytic solution analogous 
to ( p.SI ). However, we can study the dynamics of this theory analytically using perturbation 
theory. Near the false vacuum we can write 

^ = + 6(1) (3.4) 

and linearize ( |3.3[ ) in 6(j) to obtain 

(aV+5+ -1)5^ = 0. (3.5) 
This equation admits a single growing mode 

6^ = aoe^^/^""'^^) (3.6) 
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with ao constant. This solutions describes the usual tachyonic instability near the false 
vacuum. 

Near the true vacuum we take 

(l) = K^ + 6()) (3.7) 

and linearize (|3.3D in 6(j) to obtain 

(aV+9+ - 1 + 2K2a'y+a+^ ^ Q ^ 

This equation belongs to the class studied in where the generatrix is f{s) 
1 + 2K°''^^\ The full solution is 

5cP = ^ane'"'\ (3.9) 

n 

where a„ are arbitrary complex numbers and 

a'V+Sn = 1 —Wn UK^lnK] . (3.10) 

2 In K L J V ; 

Here Wn are the branches of the Lambert-W function and n runs over all integer values 
(both positive and negative). It is easy to verify that the s„ are complex and appear 
in complex conjugate pairs so that one may choose o„ to obtain a real-valued solution. 
Moreover, it can be shown that Re(s„) < for all n so that all of the modes near 
are decaying. These decaying oscillatory modes are ghost-like and (perhaps) are related to 
the presence of closed string excitations near the perturbative vacuum. 

We will see shortly that, quite surprisingly, this naively linearized analysis actually 
gives a good qualitative picture of the fully nonlinear tachyon dynamics: generic solu- 
tions roll away from the unstable maximum and undergo damped oscillations about the 
minimum. 

The reader may find it puzzling that perturbation theory around different critical 
points of the potential yields different numbers of initial conditions. This kind of mis- 
match (which is not unique to equation (^]^)) is an artifact of the perturbation theory 
employed. We discuss the resolution of this mis- match in the appendix. Note that because 
of this mis-match a naive perturbative analysis may lead to misleading results concerning 
the counting of initial data in higher derivative theories. 



(3.8) 

= aV+s - 



4. Numerical Methods 

The naive perturbative analysis employed in subsection |3.2| is, of course, not sufficient to 
establish the stability of generic solutions of equation ( |3.3D . Since we are unable to obtain 
nonperturbative analytical solutions of ( |3.3D we must turn to numerical analysis. In this 
section we describe our numerical methods. Although our primary interest is in equation 



|3.3| ), we apply our methods also to the p-adic string equation {2A) as a consistency check 



Our approach follows closely the formalism developed in [36| to study nonlocal cosmological 
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models. Ref. [36| improved significantly on previous efforts to solve nonlocal equations 
numerically by allowing the equations of motion to be solved as an initial value problem. 
Since the stability of the theory is intimately tied to the structure of the initial value 
problem, it is only in the context of this formulation that one can sensibly address the 
crucial issue of stability. 

4.1 Partial Differential Equation Formulation 

Infinite order differential equations such as (|2.4D and (3.3) are not directly amenable to 



standard numerical analysis. In order to solve these equations on a computer it is con- 
venient to introduce a fictitious auxiliary direction (which we call r) and re-formulate 
the nonlocal ordinary differential equations (ODEs) as local partial differential equations 
(PDEs) in the space spanned by the coordinates x"*" and r. (This is very much analogous 
to the "diffusion equation" formulation that has been used previously to study nonlocal 
cosmologies numerically |28, 30, ^.) 

We start by remarking that both the equations of motion ( |2.4| ) and ( |3.3| ) can be 
expressed as 



(l + 4^2^a+)e-"^+t/;(x 



where we are using the notation of Ref. 1 3£ ] . 

For the p-adic case, = 0, a = —a'V^lnp and (p = exp[—a'V^ lnp/2{p 
for the SET case we have = -1/(8 In /sT), a = 2a'V+ In/C, p = 2 and = K^e 

We now introduce an auxiliary variable r and define a new field such that 



(4.1) 



1)]^/;, and 



-rad4 



(4.2) 



By differentiating ^'(x^,r) with respect to r we find that ^{x^,r) satisfies the PDE 

1 5^'(t,r) 



a 



dr 



(4.3) 



with the boundary condition 



^'(x+,1) -4^2 



(?^'(x+,r) 
dr 



(4.4) 



which is determined from the equation of motion (44) employing the PDE (^ 

The nonlocal system has now been formulated in a form amenable to standard numer- 
ical methods. Once we have specified the initial data ^'(xj^,r) (see the next subsection) 



we can proceed to numerically integrate the PDE (4.3) subject to the boundary condition 
( |4.4| ) on the interval < r < 1, > Xj^. At the end of the calculation the solution tp^x'^) 



of the original nonlocal ODE (4.1) is extracted as 



ip^x'^ 



^'(x+,0) . 



(4.5) 



Note that the system (^^,4.4) is very similar to the diffusion- like system obtained in 
136]. In fact, it is identical with t replaced by x^, □ replaced by 9+. The diffusion-like 
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system obtained in ||36[] was ill-posed in the sense that high frequency initial data grew 
faster than lower frequency data, and hence numerical errors (which can be thought of 
as very high frequency noise) rapidly grew to swamp the real solution. Dealing with this 
numerical problem was a major part of the work in |36|. The ill-posedness arises because 
the PDE that was solved in |36| is second order in time and first order in the auxiliary 
direction. In the case at hand, however, the PDE ( |4.3| ) is only first order in light-cone time 
and as a result solving this PDE is the manner described is a well posed problem. This 
means that our numerical method is very stable, and the solutions produced are highly 
robust. 

4.2 Constructing Suitable Initial Data 

It remains to specify the initial data ^'(Xj^, r) which must be consistent with the boundary 



condition (4^). In general, choosing such initial data is nontrivial. In order to proceed 
we construct suitable initial data "^{x^,r) perturbatively for ip{x^) close to some value 
A. There is no loss of generality because our formalism allows us to fix A arbitrarily. 
Once we have approximately determined 'i?{x^,r) this initial configuration is then evolved 
numerically into the nonlinear regime. Hence our solutions are indeed fully nonlinear. 
For initial field value ip{xf') close to some (arbitrary) constant A we define 



V'(x+) = A-h(5V'(x+) . 
Linearizing (|4.1| ) in 6ip{x~^) leads to 

(1 + 4('^ad+) e-"^+(5V = + AP-A. 



(4.6) 



(4.7) 



Now ^(x^ 



A + (5^'(x+, r) where (5^'(3;+, r) = e-™^+(5^/;(x+) satisfies the PDE 

1 d 
= —5*. 

a or 



To linear order in 5ip the boundary condition (|4.4| ) becomes 

"(9(5^'(x+,r 



(5^'(x+,l) -4^' 



dr 



pAP~^6'^{x+,0)+AP -A. 



(4. 



(4.9) 



We can now solve for 5^{x^ ,r) by separation of variables. We take an ansatz of the 



form 



6^{x~^,r) = 5'4){x'^)g{r) + h{r) . 



(4.10) 



Substituting (|4.10D back into (^^) we find that the functions dependent on the auxiliary 
variable r are of the form g = e~°^^ ^ and h = h{e°'^ — l)/w^, while dip satisfies the local 
equation 



Here a;^ is any solution of the characteristic equation 



(4.11) 



(4.12) 
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and b is given by 



AP-A 



= — ; — ■ (4.13) 

In general the characteristic equation (4.12) will have many roots u;^. For each root 
of the characteristic equation we can obtain a solution 6'ipn{x~^) of equation ( [4.11| ) which, 
physically, can be thought of as a particle-like excitation near ip = A. The spectrum of roots 

corresponds to the spectrum of masses for these physical excitations. For each state 
5ilJn{x'^) there also exists a particular solution of the PDF: 5'^n{x'^-,T) = e°'^"^5^nix^) + 
bn(e"'^"'~ — l)/w^. To construct general solutions, of course, we must superpose these modes, 
and we are lead to the general solution 



(5^(x+,r) = 



5V'n(x+)e"<^ + 



(4.14) 



where it can be easily verified that the constants bn must now satisfy the relation 

AP-A 



Ebn 
7^ 



ojI pAP-^ - 1 



(4.15) 



which replaces ( fl.l3 ) when we superpose more than one mode, but are otherwise arbitrary. 
When .^^ = 0, such as in the p-adic case, the roots are given by 



aujl = \n{pAP-^)±i2TTn, 



and when ^ 0, like in the SFT case, we obtain 



1 

4^ 



when A^ Q and 



OLUJ 



1 

4^2 



(4.16) 



(4.17) 



(4.18) 



when ^ = 0. Wn again represents the branches of the Lambert-W function. For both the 
p-adic and SFT cases, if A is sufficiently close to the unstable maximum then there exists 
a real-mass mode with aw^ > 0. This mode corresponds to the usual tachyon and reflects 
the instability of the false vacuum. (There is also a decaying mode with real-mass which 
drops out of the spectrum very quickly.) In addition to this tachyonic state we also have 
an infinite tower of states with complex mass-squared. The kinetic energy associated with 
these states has indefinite sign |T3[| and hence these states are ghost-like. More precisely, 
each complex-mass state behaves as an admixture of a ghost field and a non-ghost field 
(hence we use "ghost-like" rather than "ghost"). In |36| such states were described as 



"quintoms" and we will occasionally use this term interchangeably with ghost-like. The 
dependence of the mass spectrum on the initial value A is illustrated in Figs. dH) and @). 
There we show the effective potential of the models under study and the regions where 
the characteristic equation has at least one real root and regions where all the roots are 
complex (which can be compared to the time-like cosmological cases in Ref. |48[|). 
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Figure 3: The p-adic string effective potential with p — 3. 



n =0, -1 n > , II, < - 1 




Figure 4: The effective potential of the SET case. The left panel shows the regions where the 
dynamics is initially described by a real field (solid line) or as a quintom (dashed lines) for n — 0, —1. 
The right panel shows the same regions for the remaining values of n. 



We can see from the preceding discussion and from Figs. (^) and (^, that the charac- 
teristic equation, for a given point A, has at most two real roots and an infinite number 
of complex roots. Hence, a particular solution 5^n{x~^ ,r) is, in general, complex- valued. 
For each complex root, its complex conjugate is also a root, hence, (5^* is also a particular 
solution to the PDF equation. Moreover, by suitably combining 6^n and 5^^, we can 
construct real-valued solutions. Having this in mind, we are now ready to specify suitable 
initial data for our simulations. We are particularly interested in the evolution which fol- 
lows when the field is initially in a region close to the maximum of the effective potential. 
It can readily be verified that an acceptable initial profile in the p-adic case, where the 
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maximum exists at ^ = 1, is 

^{x+) = 1 - e 5^ a„e"""?^^ cos(ac^„?r + 9^) , (4.19) 

n=0 

and an initial profile for the SFT case, where the maximum is at -0 = 0, can be written as 

^(^1^) = e Yl cos(aw„?r + 0^) , (4.20) 

n=0 

where 9n is an arbitrary phase, the indices R and / mean the real and imaginary parts, 
and e is a small number. 

In equations ( 4.19| ), ( |4.20 ) the free coefficients Un allow us to fix an infinite number 



of initial conditions. Physically each a„ parameterizes the amount of the state with mass- 
squared that is present in the initial admixture. The approximate initial functions 
( |4.19 ), ( |4.20 ) can now be numerically evolved using (|]^), ( [4.4] ) into the fully nonlinear 



regime. 

For simplicity, we have set the arbitrary phases 9n to zero in our examples. We 
have verified that the inclusion of ^„ 7^ does not qualitatively change the behaviour 
of our solutions. In particular, the inclusion of these phases does not change our results 
concerning the nonlinear stability of generic SFT solutions. 

4.3 Consistency Check Using the p-adic Theory 

Although our motivation for turning to numerical analysis was equation ( |3.3D , we can 
also use this method as a consistency check on our previous results for the p-adic string 
(section ^). First, let us verify that we can reproduce the analytic solution ( ^.5| ) using the 
PDF formulation. To this end, we introduce another field T(x^,r) related to ^{x~^,r) by 



T = In^*. In terms of this new field the boundary condition (4^) is linear: T(x^,l) = 
pT{x^,0). Now we can solve the diffusion-like equation ([4.3| ) exactly by separation of 
variables, T = f{x~^)g{r) where f{x~^) and g{r) are of the form 



/(x+) = e--'^+ , (4.21) 
,2 



<7(r) = e"-^ (4.22) 

The roots of the characteristic equation auj'^ are determined using the boundary conditions 
to be: 

auj'^ = In p ± i27Tn . (4.23) 

Putting all these results together and using ( [4.5| ) leads to the same solution for (p(x~^) 
found in Eq. (2.5). This confirms that the PDF method is consistent with other methods 



of solving non-local equations of motion. 

Note that for the SFT case, the boundary condition ( [4.4| ) cannot be written in terms of 
a linear relation, and therefore the method of separation of variables cannot be employed 
to obtain a solution to the PDF. This is simply a reflection of the statement we made 
previously that we cannot find an easy analytic solution to Eq. ( |3.3[) (unlike the p-adic 
case Eq. (^)). 
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light-like p-adic string 




Figure 5: Evolution of tlie combination of the real field with the first quintom (n = 1) in the 
p-adic string (solid line). The dashed and dash-dotted lines represent the evolution of the real field 
and of the quintom, respectively, when they are rolling in isolation. We have used e = 10"'^, oq = 1, 
ai = 0.25, p = 3, and = 0. The x+ axis is measured in units of a'V^. 



As a final consistency check we compare numerical solutions of (^^) and ( [4 .41 ) using 
the approach of subsection 4.2 to the exact analytic solutions obtained in section In 
Fig. 1^ we plot the light-cone time evolution of the p-adic tachyon for a combination of the 
real field (n = 0) with a oscillatory mode, which we refer to as a quintom in line with p6| . 
We choose the n = 1 mode. We see that the real field helps the quintom to decay and 
to remain harmless for the entire length of the subsequent evolution, assuming that the 
amount of quintom present initially does not violate the condition ( |2.8D which we derived 
analytically above. Moreover, the numerical solution is in excellent agreement with the 
analytical solution ( |2.5D , a further verification of both the interesting behaviour of that 
solution, and of our numerical method. 



5. Nonlinear Stability of the Light-Like SFT Tachyon 
5.1 Numerical Solutions 

Having described in detail our numerical methods in section ^ we now wish to apply these 
to study the nonlinear dynamics of equation light-like tachyon condensation in SFT for 
arbitrary initial data. We proceed by numerically integrating the PDF (|4.3| ) using the non- 
linear boundary conditions (4^), from initial data of the form ( [4.19D -( ^.20[) determined by 



perturbing about the hill-top as we have just described above. We stress that, although the 
initial data is fixed by considering a linearization of the equations of motion, the numerical 
solution rapidly evolves out of this linear regime and behaves in a fully non-linear manner. 

First, we verify that picking initial conditions such that the real field tachyon is present 
initially, leads to the well behaved rolling solution which was constructed in |14]. Figure ^ 
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Figure 6: Evolution of the real field n = in the SFT case. The upper panel shows the early 
evolution where we can see that the field rolls straight to the minimum of the potential. The lower 
panel shows the full evolution. 

presents the rolling solution produced by our numerical method which follows from picking 
only the n = initial condition. As can be seen, this solution has exactly the behaviour 
expected. 

We now focus on the fate of the ghost-like excitations. In Fig. |^ we show that the 
quintoms decay on their own in the SFT case without need of the extra contribution from 
the real field. Starting the evolution with only the n = 1 quintom present, we see that 
this excitation first decays towards the unstable maximum, and then at late times a real 
field is formed which evolves towards the minimum of the effective potential performing 
damped oscillations (assuming, of course, that the real field does not start to roll down the 
unbounded side of the potential). The initial decay is in line with the expectations drawn 
from a perturbative analysis close to the hill-top, then the subsequent roll to the minimum 
is the natural consequence of the maximum being an unstable point which the real field 
tachyon naturally tries to roll away from. 

We have explored a wide range of other initial conditions for the SFT case, and other 
choices of n, and the situation is always the same: any ghost-like excitation present initially 
decays, as we expected from a perturbative analysis, and at late times the rolling solution 
is reached and the field decays to the minimum of the potential. Indeed even if we fix our 
initial conditions near to the minimum of the potential by perturbing about a point near 
to the minimum (where the roots to the characteristic equation are all complex and we can 
only have ghost-like states present as an initial condition), the behaviour that is always 
found is a rapid decay to the minimum. 

5.2 Ghosts, Stability and Friction 

The dynamics of equation ( |3.3[ ) are quite remarkable. This equation admits infinitely 
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Figure 7: Evolution of the quintom n = 1 in the SFT case. The upper panel shows the early 
evolution where we can see the quintom decaying towards the maximum of the potential. The lower 
panel shows the full evolution where it is clear that only a real field evolving towards the minimum 
subsists. 



many initial conditions and, by Ostrogradksi's theorem, the Hamiltonian is unbounded 
from below. However, the dilaton gradient violates time translation invariance so that 
intuition with Hamiltonian dynamics gives a completely inaccurate picture of the actual 
behaviour of the tachyon. Before getting too excited, we should remember that this stability 
relies crucially on the ansatz (/> = (j)(x~^) and that on small scales (compared to the dilaton 
gradient) generic field profiles should display the usual unstable behaviour. With this 
caveat in mind, however, we believe that our findings our quite significant. Equation 
provides an explicit example of an interacting nonlocal theory that admits infinitely many 
initial conditions but is completely stable. It is our hope that this toy example will provide 
hints into how to construct more realistic stable theories with infinitely many derivatives. 

It is worth noting that our discussion of stability here refers only to the kinds of 
erratic time evolution that are associated with the Ostrogradski instability. Of course, the 



potential for the theory (3.3) is unbounded from below and hence there is an instability 
associated with rolling down that unbounded direction. This instability would be present 
also in a local field theory with the same potential and is completely unrelated to the 
kinds of higher derivative instabilities that we are interested in. The unboundedness of the 
potential is thought to be physically associated with the closed string tachyon. 

6. Stability of Nonlocal Cosmologies 

Owing to the close similarity between the equations describing light-like tachyon condensa- 
tion and the equations describing nonlocal cosmologies (such as p-adic inflation) one may 
wonder whether the kinds of phenomena discussed above occur also in the latter case. It 
has been observed previously that it is possible to obtain cosmological solutions where the 
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tachyon settles down to the true vacuum at late times |30, However, previous studies 



have not considered how generic such solutions are, and indeed other wildly oscillatory 
solutions have also been produced |36]. 

We have numerically studied both p-adic and SFT cosmologies for a variety of initial 
conditions describing a mixture of tachyon and ghost-like excitation in the initial state. 
Our preliminary results suggest that, at least for certain parameter choices, the set of 
initial conditions leading to non-pathological evolution is not measure zero. Hence, we 
are lead to suspect that the island of stability is a fairly general property of infinite order 
theories in the presence of friction. This suggests that the combination of friction and 
constraints on the initial data provides a very general recipe for obtaining stable solutions 
in higher derivative theories.^ This observation should be helpful to guide future attempts 
to construct physically sensible time-dependent solutions in SFT. 

Before we leave this section a few comments are in order. As discussed previously, in 
the cosmological context our numerical methods are much less robust than in the light-like 
case. Although our preliminary efforts suggest that the island of stability exists in the 
cosmological context, it is difficult to make conclusive statements because our finite nu- 
merical accuracy prevents us from following the evolution to arbitrarily late times. Hence, 
although our solutions appear to settle down to the minimum of the potential, we cannot 
(yet) rule out the possibility that instabilities re-appear at very late times. (Note that this 
problem does not afflict our light-like solutions discussed previously.) We intend to return 
to this issue in future work. 

7. Are Wild Oscillations Necessarily Catastrophic? 

Throughout this paper we have discussed the nonlinear dynamics of the infinite order equa- 
tions that describe tachyon dynamics in string theory. Our focus has been on constructing 
solutions that are not afflicted by wild oscillatory (or otherwise unstable) behaviour at late 
times. However, we have not considered the question of whether such instabilities are nec- 



essarily catastrophic. If a field theory such as (2A.) or ( |3.1[ ) were the complete picture the 
answer might be straightforward, however, the problem is rather more subtle in the context 
of string field theory. The stress tensor at a point is not a truly gauge invariant object 
since neither the interactions of off-shell close string vertex operators, nor the restriction of 
the boundary state zero mode to particular values, preserves BRST invariance |14]. Since 



it is not trivial to relate the tachyon (/> in the level truncation to physical observables it 
is not entirely clear if the wild oscillations are problematic. Several different explanations 
have been proposed in the SFT literature. We briefly discuss these below. 

Erler and Gross argued, using a clever choice of basis (the light-cone basis), that the 
full SFT is first order in a single null direction (say x~^) and nonlocal in all the remaining 



directions [49|. The initial value formulation in this case is somewhat more complicated 
than in nondegenerate second order systems (where one specifies the coordinates and ve- 
locities at t = 0). In the Erler and Gross formulation one must provide information about 



^Note that these kinds of constraints could be elegantly implemented at the level of the formal pseudo- 
differential operator by using the contour deformation prescription advocated in 
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the field at = and also on the surface = c for all x'^ (here x~ is the other null 
direction orthogonal to x^). In it was argued that one can take c — > — oo and demand 
that the field vanish in this limit. Although this requirement is physically sensible, it is not 
clear if it represents an undue restriction on the free coefficients of the solution. Note also 
that practical computations in the light-cone basis are complicated by the re-appearance 
of spurious negative energy states as artifacts of the level truncation pO|]. 

Coletti et al. |51| argued that the wild oscillations of the tachyon field can be eliminated 
by a nonlocal transformation of the form (j){t) = f{dt)T{t) + ■ ■ ■ which takes the cubic 
open string field theory action to the analogous boundary string field theory action. The 
generatrix f[s) in this case has both poles and zeroes in the complex s-plane and hence 
one might worry about changing the number of degrees of freedom in the solution (see, for 
example, |12|). 



Finally, Kiermaier et al. [52| constructed a class of BRST- invariant closed string states 
for any classical solution of open string field theory. This state can be used to provide gauge- 
invariant observables. Using this state Kiermaier et al. argue that the wildly oscillatory 
rolling tachyon solution actually describes the regular close string physics studied in Q. 
The peculiar time evolution is interpreted as resulting because the regular physics of the 
closed string sector is being described in terms of open string degrees of freedom. 

It is sometimes argued that the ghost-like modes present in level-truncated SFT and 
p-adic string theory are artefacts. In this case one expects that more realistic dynamics 
will be obtained by projecting these states out, presumably through some prescription for 
choosing initial data (see [12| for an elegant implementation). Our analysis of the island 
of stability can be thought of as elucidating the minimal constraint on the initial data 
necessary to obtain sensible evolution. The idea of projecting out ghost excitations from 
higher derivative theories using some boundary conditions is not new. Qualitatively similar 



prescriptions have been employed in the finite derivative case; see m] and [54|-|5q]. 



8. Conclusions 



Using a combination of analytical and numerical methods we have investigated the nonper- 
turbative stability of light-like rolling tachyon solutions in the presence of a linear dilaton 
background. We have uncovered some potentially surprising results. We have seen that 
the addition of friction can drastically soften the effects of higher derivative instabilities. 
In the case of the p-adic string (and also VSFT) we have found an island of stability in 
initial condition space. For initial conditions within this island the tachyon dynamics are 
non-pathological. Interestingly, the island of stability is not a set of measure zero nor is 
it an artifact of perturbation theory. We have found qualitatively similar behaviour in 
the cosmological context and have speculated that the recipe of mixing friction with some 
constraints on the initial data provides a general prescription for constructing sensible 
(particular) solutions in nonlocal theories. 

In the case of SFT at level zero truncation the effect of the friction on the higher 
derivative instability is even more dramatic. In this case the unstable growth associated 
with the ghost-like modes is completely damped out by the dilaton gradient and the re- 
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suiting tachyon dynamics are non-pathological for generic choice of initial conditions! This 
provides an invaluable example of an interacting nonlocal theory derived from string theory 
that is completely stable. A caveat is that this stability relies on the ansatz (j) = (j){x~^) and, 
as we have argued, is not expected to persist in the case where 4> depends on all space-time 
coordinates. We do not believe that this caveat should be viewed as a serious limitation 
on our analysis since the significance of our results does not lie in the claim that light-like 
profiles are the most realistic rolling tachyon solutions. Rather, our results are important 
because we have uncovered a previously unexpected loop-hole in Ostrogradski's theorem. 

During the course of our investigation we have developed many general techniques 
for studying nonlocal theories at the fully nonlinear level. It is our hope that these re- 
sults/techniques will lead to the discovery of more realistic examples of stable infinite 
order theories. 
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APPENDIX: The Perturbative Mis-Match 



In subsection 3.2 we discussed the perturbative construction of solutions of equation ( p.3| ) 
about the constant solutions = 0, . There we were presented with the curious puzzle 
that - taking this result seriously - one would infer different numbers of initial conditions 
about these two critical points. This mis-match is an artifact of the perturbation theory 
that was employed. To see this, we reconsider solving equation ( ^.2| ) in a linearized expan- 
sion about an arbitrary value (p = We could simply extract the result from our analysis 



in subsection [4.2| . However, the formlism does not rely on the PDE formulation and it 
may be of interest to show how to perturb about an arbitrary without introducing the 
auxiliary direction r. 

We begin in a very general context and specialize to equation (^]^) at the end. Consider 
the nonlocal equation 

= V [cl>] (A-1) 

where F{z) is an analytic function of the complex variable z that can be represented by a 
convergent series expansion 

oo 

F(z) = ^a„z" (A-2) 

n=0 

and D is some linear differential operator which satisfies D(const) = 0. We wish to solve 



( A-1 ) near some constant value (p = (j)Q not necessarily a solution of ( A-1 ). Writing 

^ = ^o + 6(j) (A-3) 
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and linearizing in 5(j) we have 

[F{D) - V" [M] Sep + [F(O)0o - V [M] = (A-4) 

If (j)Q were a solution of ( \A-l\ j then the second term in the square braces would vanish. In this 
case one could construct the solution 5(j) by provisionally taking S(f) to be an eigenfunction 
of D (that is assuming Ddcp = —uj'^dcp). This is precisely the approach that was adopted 



in |^^-[^4[ and |37|. However, we can easily generalize this approach more general values 
of 00 by provisionally taking Sep to be a solution of the equation 

D6(l) = -uj'^6()) - b (A-5) 

It is straightforward to show that the action of the pseudo-differential operator F{D) on a 
function satisfying ( [A-5| ) is 

F{D)6(I) = F{-oj^)6(l) + [F{-uj^) - F(0)] (A-6) 

Plugging ( |A-6D into ( |A-4| ) we find that the solutions of ( |A-5D are also solutions of the 
fully nonlocal equation ( |A-4D as long as 6 are chosen to satisfy the following algebraic 
equations: 

= F(-u,^-)-V"\<t,o] (A-7) 
"' vm-Fm (A-8) 

In general equation ( [A-T]) may have multiple roots u>n {n = I;""" lA^) which defines 
solutions bn of equation ( A-8| ). There are then N independent solutions dcpn (each solving 
an equation of the form D6(f)n ~ ^n^'t'n bn ) and the most general solution of ( [A-4| ) is 
obtained by superposing these modes as 

AT 



J2 ^'^n (A-9) 



n=l 



When superposing multiple roots equation ( |A-8D generalizes to 

^bn _ -F(O)<^o + r[0o] 

^ V" [M - F{0) ^ ' 

Note that this method may fail if = is a solution of ( |A-7 ). This approach is identical 
to the formalism employed in and also in subsection |4.2| . 

Now we specialize to equation (p^). We take D = 9+, F{z) = [a'V^z — i^^^-Sa'v^+z 
and V' [(/>] = —(fP'/K^. After some straightforward manipulations we find that the mode 
functions (50„ take the form 

50n(x+) = a^e^"^'"- + ^ (A-11) 



where the are arbitrary constants and 

2\nK K 



(A-12) 
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(where Wn denotes the branches of the Lambert-W function and n runs over all integer 
values). Summing over the modes and making use of equation ( A-IOD we have 



n 

For = these expressions reproduce equations ( p. 91 ) and ( P-lOp . Let us consider the 
case (/>o = 0, where a naive perturbative analysis yields only a single growing mode. Taking 



the limit 0o ^ in equation ( A-12 ) we find that sq l/{a'V~^) while Re(s„) — > — oo 
for all n 7^ 0. Hence, the infinite tower of ghost-like modes decay very quickly near the 
false vacuum. In the limit (j)Q ^ these spurious states all go to zero infinitely fast and 
completely drop out of the spectrum. Thus, the analysis in this appendix is consistent with 
equation (|3.6| ) and the puzzle of the mis-match of initial data counting is resolved. 

A similar mis-match would occur if we have studied the p-adic equation ( |2.4| ) in a 
perturbative expansion about the constant solutions (j) = 0,p~"^ /[2{p-i)] There one 
finds infinitely many solutions about the false vacuum and no non-trivial solutions at the 
true vacuum. Again, the mis-match is an artifact as can be seen by examining the fully 
nonperturbative solution (|2.5| ). From this solution it is clear that the extra states become 
nonperturbative near (j) = 0. 
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